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Abstract 

Radion stabilization is analyzed in 5-dimensional models with branes in 
the presence of Gauss-Bonnet interactions. The Goldberger-Wise mecha- 
nism is considered for static and inflating backgrounds. The necessary and 
sufficient conditions for stability are given for the static case. The influence 
of the Gauss-Bonnet term on the radion mass and the inter-brane distance 
is analyzed and illustrated by numerical examples. The interplay between 
the radion stabilization and the cosmological constant problem is discussed. 



1 Introduction 



The idea that we hve on a brane in a higher dimensional space-time has 
attracted a lot of attention since Hofava and Witten presented [T] their 
1 1-dimensional model with 10-dimensional branes, motivated by M-theory. 
Many 5-dimensional models with 4-dimensional branes (some of them being 
effective models derived from the Hofava-Witten one) have been discussed 
since then. In any such model, the distance between the branes should be 
stabilized. A scalar field related to that distance is usually called the radion. 
Several, more or less precise, definitions of the radion are used in the liter- 
ature. Of course, not every such a definition is reasonable. It was pointed 
out j2] that the radion should at least fulfill appropriate linearized Einstein 
equations of motion. It is easy to check that in models where only gravity 
propagates in the bulk the radion corresponds to the only gauge-invariant 
scalar perturbation of the metric. Moreover, such a radion is massless, which 
means that the distance between branes is not stabilized. Goldberger and 
Wise |H1 proposed a class of simple models to stabihze this distance. They 
introduced a bulk scalar field with interactions described by some bulk and 
localized brane potentials. Perturbations of such a field mix with the scalar 
perturbations of the metric, giving rise to an infinite tower of KK states in 
an effective 4-dimensional theory. The radion can be defined as the lightest 
of those states. This lowest KK mode might sometimes have negative mass 
squared, which means that a given configuration is unstable. Situations 
where such a tachyonic radion may appear were discussed in several papers 
in El El- In quite general criteria for such instabilities were presented. 

One of the possible generalizations of this type of models takes into ac- 
count interactions of higher order in the curvature tensor. Such interactions 
naturally appear in the a' expansion in string theories [7j. The simplest 
correction is the famous Gauss-Bonnet (GB) term. Many aspects of 5- 
dimensional models with GB interactions have been already discussed in the 
literature. However, the influence of GB interactions on the radion stabiliza- 
tion has not been addressed^ . It seems that an appropriate way of modeling 
the radion is by introducing an additional bulk scalar field a la Goldberger 
and Wise. In this work we analyze the radion stabilization in models with 
GB interactions and with such a bulk scalar field. We consider 5-dimensional 
backgrounds with flat static 4-dimensional sections as well as with inflating 
ones. The case with inflation is especially interesting, because it was shown 

^The problem of radion stabilization in models with the GB term has been recently 
addressed in |H]. However, the definition of the radion used in that paper does not seem 
to be adequate. It is considered there as the (5,5) component of the metric. A potential 
for such a field is introduced and this explicitly breaks symmetries of general relativity. 
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P that during inflation stabilization of the radion is more difficult. We find 
some generic conditions for radion stability in terms of background solutions. 
We discuss these conditions using analytical arguments as well as results of 
numerical calculations. 

Models with scalar fields and Gauss-Bonnet interactions have been con- 
sidered in the literature, but in contexts other than radion stabilization. 
Problems of singularities and fine tuning in such models were discussed in 
|inj . Phantom cosmology with GB corrections was investigated in [TT| I12j. 
In jini, the interplay between the GB term and the quintessence scalar was 
analyzed. Experimental constraints on the quintessence- GB coupling were 
derived in Cosmological models with scalar dependent GB interactions 
were considered in [T3]. In jTH], a scalar field was used to support a smooth 
brane in an Einstein-GB gravity model. A cosmological model with GB 
interactions and two scalar fields was investigated in [T7j . 

The organization of this paper is as follows: In the next section, we present 
the equations of motion and the boundary conditions for models with GB 
interactions and the Goldberger-Wise scalar field. The results are presented 
for an ansatz describing the 5-dimensional space-time with 4-dimensional 
inflating de Sitter sections and two branes. In section 3, equations of motion 
and boundary conditions for the scalar perturbations in such a model are 
derived and discussed. A variational formula for the radion mass is given in 
section 4. It is used later to derive the necessary and sufficient conditions for 
radion stability in the case of a vanishing Hubble constant. In section 5, we 
compare two classes of solutions present in models with the GB term. It is 
shown that the stability analysis strongly disfavors the so called "new" so- 
lutions, which appear to be always unstable. Section 6 contains a discussion 
of the changes caused by the presence of GB interactions in the action. The 
relation between the radion stabilization and the cosmological constant prob- 
lem is shortly discussed. Some results of numerical calculations are presented 
in section 7. Section 8 contains our conclusions. 



2 Background solutions 

We consider 5-dimensional models described by the action 



S 



r 2 



where is the Gauss-Bonnet term quadratic in the curvature tensor: 

- AR^.R^'' + R^,p„R^'''\ (2) 



r2 



'^fj.upcr-' 
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The scalar field $ interacts via the bulk potential V{^) and two brane po- 
tentials ?7j($) located at yi and 1/2. We do not include extrinsic curvature 
terms in the action, because we work in the up-stairs picture of the 3^/7^2 
orbifold, where all total derivatives integrate to zero. The boundary condi- 
tions are obtained by integrating the equations of motion around the branes. 
This method is equivalent to using the Israel junction conditions |18jj after 
adding appropriate extrinsic curvature terms to the action [T^ . 

We are interested in warped background solutions with the following 
ansatz for the metric and the scalar field: 



ds2 = a{yY {-dt^ + e^^'6ijdx'dx^ + dy^) 
$ = (P{y). 



(3) 
(4) 



It describes the 5-dimensional space-time warped along the ^/-coordinate 
with the 4-dimensional de Sitter slices inflating with the Hubble constant 
equal to H. For H = 0, the same ansatz describes a warped geometry with 
fiat Minkowski foliation. For simplicity, the 3-dimensional space has been 
chosen to be fiat. Also for simplicity, we assume that inflation is driven by 
some brane fields which do not directly influence the bulk scalar perturba- 
tions, and that the inflaton and radion sectors can be treated separately. 
In this paper, the Hubble constant is treated just as a given parameter. Of 
course, analysis going beyond such a simplification would be very interesting. 

The bulk equations of motion for the system described by action and 
satisfying ansatz (j^HH) are given by (we use units k = 1) 



f + 3-0' 
a 



aV = , 



a \ a 



A'2 



0, 



1 



1 + ^1- ^0'^ + a'V = 



(5) 
(6) 

(7) 



The primes denote differentiations with respect to the appropriate (implicit) 
arguments, i.e., with respect to the scalar field $ in the case of V, and 
with respect to the coordinate y for a and (p. To simplify the formulae, a 
new function - depending on the warp factor, the Hubble constant and the 
coefficient of the GB term - has been introduced, namely 



^ = a -4a 
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The modifications of the background equations of motion caused by the GB 
term can be writen in a simple way in terms of just this single function. It 
is interesting that, as we will show later, this also holds for almost all other 
modifications (equations of motion for the perturbations, stability conditions, 
most of the boundary conditions, etc.) considered in this paper. The only 
exception is one of the boundary conditions we discuss below. 

The boundary conditions are obtained by integration of the equations 
of motion over infinitesimally small intervals around the branes. Of course, 
one must add to the eqs. (jSHZj) the contributions coming from the brane 
potentials Ui. Such integration of eq. © gives the boundary condition for 
the derivative of the scalar background: 



lim 



2 ' 



(9) 



The corresponding condition for the warp factor a{y) can be obtained from 
eq. ^ after using eq. ((Tj) to eliminate 0'. The result reads 





o 



(10) 



The two above equations describe the jumps of the Z2 odd functions a'{y) 
and (f)'{y) at the Z2 fixed points yi and ?/2- Observe that the expression in 
the square bracket of the last equation differs from ^ defined in (jHl) just by a 
factor 1/3 in front of the term containing a'^. The origin of this additional 
factor is the following: As the derivative a' is a function discontinuous at 
the positions of the branes, the second derivative a" contains Dirac delta 
contributions. Hence expressions of the type a"{a'Y must be regularized 
before integration. Any regularization leads to the same result 



'^^\"{a'fdy=\{a'f 



(11: 



which gives the factor 1/3 in eq. |TO|l . 



3 Scalar perturbations 

The main interest in our analysis is the dynamics of scalar perturbations 
around the background metric given by the solutions to the equations of 
motion (jSHZj) with the boundary conditions (jn HTIH) . In general, one could 
perturb all the components of the metric, the scalar field and the positions 
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of the branes. Some of such perturbations correspond to gauge degrees of 
freedom. This gauge dependence has been carefully discussed in the literature 
1201 EH |22j . One can choose gauge independent combinations of perturba- 
tions or work in a specific gauge. In the following we will use the generalized 
longitudinal gauge. The scalar perturbations can be written in this gauge as 

ds^ = [(1 + 2Fi) (-dt^ + e^^'dijdx'dx^) + (1 + 2F2) dy^] , (12) 



$ = + ^3, (13) 

where a and (p are the y-dependent background solutions to eqs. (jSHZ|) and 
(pl HTTH) . The perturbations Fi are functions of all the 5 coordinates. This 
new ansatz has to be substituted into the equations of motion obtained from 
action and expanded in powers of the (small) perturbations Fj. The 
zeroth order equations are of course given by (jHHZI) and (j9 HT0|) . The equations 
linear in the perturbations can be used to determine the mass eigenstates 
for such perturbations in a given background. Higher order terms of the 
expansion correspond to interactions which are more than bilinear in the 
perturbations. We are mainly interested in the masses of the 4-dimensional 
scalars, especially the radion, so we only need the linearized equations. 

The off-diagonal components of the linearized Einstein equations give 
the following conditions, which must be satisfied in order to stay in the 
longitudinal gauge: 

+ -F2 = 0, (14) 

(eFi)' + ^aVF3 = 0. (15) 

The diagonal Einstein equations, combined with the background equations 
of motion (jSHZj), give the dynamical equation for the scalar perturbations: 



1 

n2 



I / /\ 2 



{U + AH^)F,+A'^F[-Ar-] F, 



3 

1 ,„ a' A „ 1 



^^.r^^-4>')F.--^4>'n = 0, (16) 

where □ is the 4-dimensional D'Alembert operator on the de Sitter slice 
of the 5-dimensional space-time (which will become the ordinary flat space 
D'Alembert operator when we later set the Hubble constant H to zero). The 
boundary conditions for the scalar perturbations read 

hm (F^ - F20') = . (17) 
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Equations (fT^ and (fT^ show that in the longitudinal gauge the three 
scalar perturbations are not independent. These equations can be used to 
eliminate F3 and F2, and to rewrite the dynamical equation of motion (|T6|) 
in terms of Fi alone. In addition, taking into account the symmetry of the 
warped background (jBHlI), we can separate the variables in the equation of 
motion for Fi. Defining 



F,{t,x,y) = Y,F^2{y) J d'kf^^2^k){t)e''' 
we obtain the following equations 

f{m?,k) + 3-ff/(m2,A:) + [a'^^^k'^ + ) /(m2 = 0, 



(19) 



-^m2 + 


4 


a' 
a 


0"! / 






+ 




(a 









F^2 = 0, (20) 



where the separation constant turns out to be the mass squared of scalars 
in the effective 4-dimensional description. 

Let us now rewrite the boundary conditions (fT7j) in terms of F^2. The 
l.h.s. of ()17|) describes limits (from below or from above) of some expressions 
smooth in the bulk. As a result we can rewrite those expressions using the 
bulk equations (fT^ and (fTH|) in order to replace F2 and F3 with Fi. However, 
this way the second derivative of Fi appears in the boundary conditions. It 
can be eliminated with help of the dynamical bulk equation ()20j) . The final 
result takes the form 

±6i(2) lim (f^^ + jFm^] + [rr? + m'') lim (F^^) = , (21) 
where 

bi = liin daU'^ + b2=\un (\aU'i - - + (^2) 

y^y+\2 a (j)' J y^y^K^ « / 

The above boundary conditions become particularly simple for a new 
variable Qm? = C-^m2 (we omit the subscript m? to simplify the notation): 

hQ'{yt) + {m^ + AH^)Q{yt) = ^, 
-b2Q'{y2) + (rn^ + ^H^)Q{y2) = 0. (23) 
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The dynamical equation of motion for Q has also a very simple form 

-{pQ')' + qQ = \vQ, (24) 



where p = 3/(2a0'2), q = {a^^')/{2a'^^) and A = (m^ + 4H^). This is the 
standard form of the Sturm-Liouville differential equation. The boundary 
conditions can also be written in the usual (see e.g. [22]) form 



where in our case cr(?/j) = — A/6j. The d/dn differentiation is in the direc- 
tion of the outer normal at the boundary. For our 1-dimensional equation 
it is (— d/dy) at yi and {+d/dy) at y2. The above boundary conditions 
are self-adjoint, but unfortunately quite unusual, because a depends on the 
eigenvalue A. This dependence is caused by the procedure of obtaining the 
boundary conditions described before eq. (PT|) . 

With such unconventional boundary conditions (j25|) depending on eigen- 
values, it is not possible to apply directly the analysis of the Sturm-Liouville 
systems from the standard textbooks [23J. However, it is possible to modify 
the standard arguments to obtain interesting results for the case at hand. 



For the problem of the inter-brane distance stabilization, the most interest- 
ing feature of the analyzed Sturm-Liouville equation is its lowest eigenvalue 
-^0 = + 4if^. There are three possibilities. First: the lowest eigenvalue 
corresponds to positive ttiq, which means that the background solution is 
stable against scalar perturbations, and the positions of the branes are really 
stabilized. Second: the lowest eigenvalue gives vanishing mg, the radion is 
massless and the inter-brane distance is not stabilized. Third: there is at 
least one tachyon signaling the instability of a given background configura- 
tion. 

We start our analysis with identifying a variational problem which is 
equivalent to the differential equation with the boundary conditions 
(j2Sl). Multiplying eq. by Q and integrating over the 5-th coordinate, 
one gets 



dQ 

dn 



(25) 



4 Stability conditions 
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where the subscript A is added to stress the relation between an eigenstate 
and its eigenvalue. This equation can be used to express the eigenvalue A in 
terms of the corresponding eigenfunction Qx as follows^: 



The smallest eigenvalue, Aq, can be obtained by minimizing the r.h.s. of the 
above equation over all smooth functions Q defined on the interval [2/1,2/2] 

Using the relation between A and one obtains the following formula 
for the radion mass 

For a general background (solution to eqs. (jHHZI)), the above expression can 
not be minimized explicitly. But we can obtain some bound using just one 
trial function Q, namely the constant one. Recalling the definitions of p and 
q, we get 

In the limit of 6j ^ 00 (very stiff brane potentials, i.e. potentials with very 
large second derivative) and ^ (no GB interactions), this expression 
simplifies to an analogous bound obtained in [HJ. 

In models with inflating branes, i.e. for > 0, stability of the inter- 
brane distance occurs when Aq > 4if^. The formula (j^Hj) can be used to 
check this condition, but in practically all models this can be done only by 
numerical calculations. However, some interesting analytic results can be 
obtained for the static branes for which A = m^. In such a case, the question 
of stability reduces to the issue of the sign of Aq. In the rest of this section, 
we assume H = 0. 

It will prove useful to rewrite the differential equation (j2l|l in terms of 
other variables. Using the definitions 



^In the case of more conventional, A-independent boundary conditions, the contribu- 
tions from the boundaries appear with the opposite sign in the numerator 
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and the background equations (jSHZj), the equation of motion 
written in the following simple form 



u"+ X 



1 



u 



0. 



can be 



(32) 



The differential equations ()2()j) . ()24j) and (jHSI) are not very convenient for 
background solutions with 0' vanishing at some y g]?/i,?/2[- The reason is 
that some of the coefficients in those equations diverge for 0' = 0. In such 
cases it is better to use a generalization of the Mukhanov variable [21] , defined 
in our model as 



,3/2 



V = a 

for which the equation of motion reads 



-^3 7-^1 



v"+{\ 



(1/^) 



0. 



(33) 



(34) 



Using the background equations (j^HZI), one can show that 



Ml! 



15a'2 



6 ^ 



1 a^f 0'^ 2 aV'* 4 a^0'\/' 

3 ^2^/ + 9 ^2^/2 + 3 ^fl/ 



2t/// 



(35) 



which is explicitly regular for vanishing 0'. The variables v and Q are related 
by 

\Q = -a^l^(t)'e~\ev)' . (36) 



Now we check if there is a massless scalar in the model. Equation ()32|1 is 
simple enough^ to be solved for vanishing A almost explicitly: 



UQ{y) = 9{y) 



CI + C2 !\y'e-\y') 



(37) 



Assuming that 0'(?/) 7^ for all y in the bulk, one can rewrite ()37p in terms 
of Q. Using also the boundary condition (j23|) at the first brane (for non- 
vanishing 61), one finds the following result 



a'{y) 



a{y) 



yi 



a'iy') 

a\y')i'{y') 



Vi 



a'{y>) 



(3^ 



(39) 



■^Unfortunately, for non-zero Hubble constant H the equations of motion for the vari- 
ables u and V can not be written in such a simple form as H32I) and H34|l . 
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where the normahzation has been chosen as Q{yi) = 1. The boundary con- 
dition at the second brane, &2<5o(2/2) = 0, can be fulfilled in several ways. 
It seems that one of the possibilities could be the vanishing of the integral 
j^^[a?^' / a'). This, however, can not be realized because the integrand func- 
tion must be positive for all y. The reason is that this function is closely 
related to the kinetic term for gravitons. This can be seen by considering 
(traceless and transverse) tensor perturbations Fj^J of the background met- 
ric 0. By expanding action to the second order in such perturbations 
we can obtain the appropriate kinetic term. We have found that, up to 
normalization and some total derivatives, this kinetic term is given by 

(VFJ/)' . (40) 

The coefficient in front of (V-Fj^-^)^ must be negative, because only then one 
can define a tower of 4-dimensional KK graviton states with the standard 
sign of the kinetic terms (an effective action for the 4-dimensional KK states 
of scalar perturbations was obtained in ^3). Hence for all y 

iM>o. (41) 

a'[y) 

because otherwise some of the graviton KK states become ghost-like. 

Let us now return to the discussion of the boundary condition for solutions 
given by eqs. (|nH|) and (jH^ . The positivity of the ratio ^' /a' means that the 
integral ^^^{a?^' / a') > 0. Thus the boundary condition at the second brane 
can be fulfilled only when 62 or 4>'{y2) vanishes. 

One can repeat the above reasoning starting from the boundary condition 
at the brane located at 1/2. Choosing the integration constants in (jHZj) in an 
appropriate way (for non-vanishing 62)5 one finds that the massless mode 
exists in such a situation only when bi or (p'{yi) is zero. Putting both cases 
together, we find that for 0' non-zero everywhere in the bulk the necessary 
and sufficient condition for existence of a massless mode is 



&i&20'(2/i)0'(y2) = 0. (42) 

Now we identify conditions sufficient for the stability of the considered 
warped space-time with branes. From expression (j?fjl it is obvious that all 
the eigenvalues are positive if p{y), q{y) and bi are positive. This conclusion is 
correct when the functions p{y) and q{y) are finite for all y, because otherwise 
our Sturm-Liouville system becomes singular, and a more careful treatment 
is necessary (quantities bi can be divergent, in the limit bi ^ 00 the boundary 
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conditions reduce to Q'ijji) = 0). The function p{y) diverges at the 
points where (p'{y) = 0. Using the definition of ^, we get 




(43) 



which shows that q{y) diverges at the points where ^{y) = 0. The finiteness 
requirement for p{y) and q{y) is equivalent to the requirement that 0' and 
^ are non-zero for all y. However, from the background equation of motion 
(jni) we see that 0' must be zero at the points where ^ is zero. Hence the 
condition p{y) ^ ±oo is stronger than the condition q{y) ^ ±oo. 
Thus the brane system is stable if, for all y, 



These are the sufficient conditions for the stability, but we show below that 
they are in fact also the necessary ones. 

One of the above conditions, namely the positivity of C,' (y) / a' (y) , we have 
already discussed. Any non-positive value of this ratio leads to some ghost- 
like states in the KK tower of the 4-dimensional gravitons. As a result, 
the second condition in eq. is anyway necessary for the stability of the 
model. 

Below we will use the following obvious principle of variational calculus. 
Let us consider two problems of minimalization of some functionals. If for 
every function the functional in the first problem is not bigger than the 
functional in the second problem then the first minimum can not be bigger 
than the second one. 

Let us consider the following two problems. Both consist in minimizing 
the expression given in (j28|l with the same positive p{y) and q{y) but with 
different 6,. Let the values of bi and/or 62 be bigger in the second problem. 
Bigger bi means smaller (at least not bigger) denominator in (j^H|) and thus 
bigger (at least not smaller) the whole expression to be minimized. Applying 
the above mentioned rule of the variational calculus, we see that the lowest 
eigenvalue in the problem with bigger bi is bigger (at least not smaller) than 
the lowest eigenvalue in the problem with smaller bi. The lowest eigenvalue, 
Ao, is a monotonic, non-decreasing function of bi. 

We apply the above reasoning to negative bi or 62- The monotonic char- 
acter of Ao as a function of bi means that the lowest eigenvalue for any 
negative bi (62) can not be bigger than the one for vanishing bi (62), which 
was previously shown to be zero. We get Aq < for 61 < or 62 < 0. On 
the other hand, Aq can not be equal to zero, because we have shown that a 



bi>0 



>0, 



<P'{y)7^o. 



(44) 



a'{y) 
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zero mode exists only if one of the bi vanishes (we consider a situation with 
(p'iUi) 7^ 0) otherwise there is a zero mode for arbitrary bt). A background 
with a negative value of any of the boundary parameters bi is unstable. 

We have thus shown that the first two conditions in eq. ()44|1 are not only 
sufficient, but also necessary for the stability. We will now prove that the 
same is true also for the third condition, which means that all backgrounds 
with (j)' vanishing somewhere in the bulk are unstable. Let us assume that 
(f)'{y) = for some yi < y < y2- One can prove that in such a case there 
is at least one negative eigenvalue in the problem defined by the differential 
equation with the boundary conditions (j^ . 

In order to show this, we consider the behavior of the solutions to eq. pi)) 
for large negative A. Such behavior can be found using the variable v defined 
in (jHHj) . The differential equation for v is regular also for vanishing 0', so the 
limit A — s> — oo can be obtained simply by dropping the ^-dependent term in 
eq. The corresponding solution reads 

v-ooiy) = csexp[\^{y - y^)] + C4exp[-\^{y - yi)] . (45) 

Using the relation ()36|). we obtain expressions for Q{y) and Q'{y) in the limit 
of large negative A. The boundary condition at yi can be fulfilled if 
C3 ~ C4 in the leading order in l/\/—X. Hence C3 is not small compared to C4, 
and away from the first brane the solution is dominated by the exponentially 
growing term 

Q-oo(2/) ~ ^^== exp[V-X{y -yi)]. (46) 

The important feature of this solution is that it is proportional to (p'ly) 
so it vanishes at y (and changes sign if (p'{y) changes sign). This should 
be compared to the behavior of Qq. As eqs. dHHl) and imply that the 
derivative Qoiy) is positive, Qo{y) is also positive for all y. The solution 
Q\{y) changes continuously with A changing between and —00. Hence 
there must exist a negative A for which has a zero point but is nowhere 
negative. It is easy to see that such a zero point must he at y = y2, and 
that Q'^{y2) < 0. (Both features result from the simple fact that, due to 
the equation of motion, Q and Q' may vanish at the same point only for 
Q vanishing everywhere.) The above arguments may be summarized in the 
following form: if 0' vanishes at any point between the branes, there is a 
negative A for which 

Qo(y2)>0, Qx(?/2)=0, QK?/2)<0. (47) 

Thus (for non- vanishing 62) neither A = nor A = A is an eigenvalue of our 
problem, because in both cases the boundary condition at the second brane 
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reduces to Q'{y2) = and is not fulfilled. It follows also that the sign of the 
l.h.s. of the second condition in (j^Hj) is different for A = and for A = A. 
Thus there must be at least one value of A between A and for which the 
l.h.s. of that equation vanishes. Such negative A is an eigenvalue of our 
problem. Hence the radion becomes tachyonic for any background solution 
for which 0' vanishes anywhere between the branes. The third condition 
in ()44|) is necessary for the stability. This completes the proof that all the 
conditions in ()44|1 are not only sufficient but also necessary for the stability 
of our brane model. 

This is a generalization of the results obtained in The generaliza- 
tion is twofold. First: we have taken into account the possibility of GB 
interactions present in the action. Second: even in the case of the standard 
Einstein-Hilbert action we were able to analyze a bigger class of possible 
backgrounds. The authors of P assumed that a'{y) ^ for all y, while we 
have been able to show that this assumption is in fact not necessary. 

Let us discuss the last point in some detail. The equations of motion for 
the scalar perturbations Q{y) are regular at the points where a'{y) = 0. The 
presence of a' in the denominator in the definition of q{y) may be misleading. 
Equation shows that q{y) = (a^^')/(2a'^^) is in fact finite for a' = (at 
such points ^ = 7^ 0). Thus for all backgrounds for which 0' and C, are 
non-vanishing for all y, the Sturm-Liouville system defined by eqs. (j^^ and 
is not singular and our discussion of the sign of its lowest eigenvalue is 
valid also if a' vanishes at some point (s). 

5 "New" versus "old" solutions 

Action (PJ contains the Gauss-Bonnet term whose strength is parameterized 
by a constant a. The results for models with Einstein-Hilbert gravity can 
be obtained from the formulae presented in this paper by setting a = or, 
equivalently, by replacing ^ with a^. However, one should be careful because 
substituting a by zero is not necessarily equivalent to the limit a ^ 0. The 
reason is quite obvious: Taking, for example, equation ((Tj), we see that it is 
linear in the combination [(a'/a)^ — H^] for a = 0, but is quadratic in this 
combination for any a ^ 0. For any non-vanishing a, there are two solutions 
to this equation (for fixed values of (p and (f)'). One of them converges to 
the a = solutions in the a ^ limit. Such type or branch of solutions 

^In [H], a different method for analyzing the equations of motion for the scalar per- 
turbations was used. Some features of the phase space associated with the equations of 
motion were analyzed to find conditions for the lowest eigenvalue. The method used in 
the present paper seems to be simpler and easier to generalize to other models. 
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is sometimes called the "old" one. There are also "new" solutions which 
diverge for a 0, but are well defined for any a 7^ 0. We will now show 
that all the solutions from this "new" branch are strongly disfavored by the 
stability requirements. 

Any background solution must fulfill three equations at each brane. These 
are the two boundary conditions © and (fTUj). as well as the bulk equation 
(|7j). Fixing the normalization of our metric to be a(?/i) = 1 and eliminating 
(p'iyi) with help of ©, we get the following two equations 

6 {[a'{y,)f - H') (l - 2a {[a'{y{)f - H^)) = - V{ct>{y,)), (48) 



l-4a Q[a'(yi)f 



4 

1 



r/i(0(yi)) (49) 





for the two quantities a'(?/i) and (p{yi). For general potentials V{(f)) and 
Ui{(f)), this set may have several solutions which can be found only numeri- 
cally. However, quite interesting results can be obtained analytically in the 
approximation of a stiff brane potential of the form 

UM = lf^ii^-Vif + ... (50) 

with very large /ii. In such an approximation, the difference 0(yi) — vi 
is small for each solution. This small quantity appears linearly in U[ and 
quadratically in Ui. Hence in the leading order in l//ii the possible values 
of a' at the brane are given approximately by the solutions of 



a\yi 



l-4a{^[a\y,)r-H' 



]^U,{v^)^0. (51) 
o 



The l.h.s. of this equation is cubic in a'{yi) and has two extrema. It is easy 
to calculate that those extrema are at the values of a'{yi) for which C,{yi), 
being a quadratic function of a'{yi), vanishes^. Thus there is at most one 
solution with positive ^(yi). This is an "old" solution which reduces to the 
a = solution in the a — limit (if it exists for the given parameters). The 
"new" solutions behave like a'{yi) ~ in the small a limit, and always 

correspond to negative iiyi). 

The above result is very important for the stability of the model. By eq. 
(jini), the ratio ^' /a' is typically negative for negative ^, which leads to an 
instability of the model. Comparing eqs. (jH)) and (jlHjl . we see that positive 

^The factor of 1/3 present in H51|l and discussed after eq. H1U|) is crucial for this relation 
between the extrema of l|ET|) and the sign of ^. 



15 



S,' /a' and negative ^ are possible only for large negative a and large H. 
Large values of a are not realistic from the phenomenological point of view. 
Moreover, equation p9|) shows that large values of the Hubble constant tend 
to destabilize the model even more. 

We have shown above that in the stiff brane potential approximation the 
"new" solutions are unstable. Relaxing this approximation is not likely to 
improve the stability of such solutions. It was argued in the previous section 
that the radion mass squared is a monotonic function of the parameters hi. 
Equations ()22|1 show that the hi decrease when we go away from the stiff 
brane potentials approximation. We conclude that all the "new" solutions 
are strongly disfavored by the stability conditions. 



6 Role of Gauss— Bonnet interactions 

In the previous section it was shown that the additional "new" solutions, 
which appear after including the GB term in the action, are probably all 
unstable. Only the "old" solutions can have a radion with a positive mass 
squared. These "old" solutions converge to the a = (no GB term) solutions 
in the a — >■ limit. This does not, however, mean that the addition of GB 
interactions is unimportant. The numerical analysis shows that solutions 
with small but non-zero a can substantially differ in some aspects from the 
solutions with a = 0. The main features of the solutions can be understood 
analytically in the stiff brane potentials approximation. For simplicity, we 
again consider the case with a vanishing Hubble constant H. 

We start with the boundary conditions at the first brane. As discussed in 
the previous section, the value of a'iyi) is given by the solution to eq. (jSH). 
For small a, it is 



1, 



6 



a 



a'{yi)---Ui{vi) 1 + —Ut{vi) . (52) 



27 



Substituting the above result into eq. ((Tj), we find 



1, 



3 



^ -^f/i (^i) l + ^U',{v^) +2V{v,). (53) 



54 



The last two equations show that the absolute values of a'{yi) and 4>'{yi) 
grow with a. The change of a and with y is faster (slower) for positive 
(negative) a compared to the a = case. One can see that the same is true 
also away from the branes. For sufficiently big y the equations of motion 
(jSHZj) are dominated by large derivatives and the contributions from the bulk 
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potential V become subdominant. In this regime, eqs. (0) and ((7j) give 



^'(y) ^ ^ ' (54) 

Expanding the r.h.s. of the last equation in a, one can see that the back- 
ground fields change faster for positive a and slower for negative a. As the 
expansion parameter is proportional to a/a^ the expansion breaks down for 
sufficiently small a. The behavior of the r.h.s. of for small a{y) depends 
crucially on the sign of a. For negative a, it tends to a constant, and a{y) 
can be arbitrarily small. For positive a, there is a minimal value of a{y) 

2 

«min ~ gttC^ , (56) 

for which a' is still real. It is easy to find the behavior of different quantities 
when a — > amin, namely: a" — oo, ^ — 0, ^'/a' — >• +oo. A solution with 
positive a ends (if there is no brane before, i.e. at a smaller value of y) at 
a singularity with a finite value of the warp factor. The solution can not be 
extended beyond such point. 

The solution for a = has a different behavior. It ends at a singularity 
for finite y at which a — > 0, a' — > const, 0' —>■ —oo. Solutions with negative 
a also end at the singularity a = 0, and the position of this singularity 
increases with increasing However, the position of the second brane can 
not be arbitrarily close to such a singularity. The reason is that the ratio 
^'/a' becomes negative before the position of a singularity, and the model is 
unstable if ^'/a' < somewhere between the branes. The behavior of C,'/cl' 
can be obtained from the following expansion 



y = 2a' (l - 4„- - -at, + -.-^V j + O(a^) . 



(57) 



For sufficiently big y, this expression is dominated by the (f)'"^ term which for 
small a grows approximately as a~^. It has the same sign as a {a' < because 
by definition the first brane is the positive tension one). Our numerical 
results confirm that indeed ^'/a' becomes negative before a singularity for 
a < solutions. For each background solution, there is a maximal possible 
distance between the branes. For a > 0, it is given by the position of a 
singularity, while for a < it is given by the requirement that a' must be 
positive. The precise value of such a maximal distance depends, of course. 
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on the details of the potentials present in action This maximal distance 
between the branes is a decreasing function of a > 0. For negative a, it is 
more model dependent. 

The positivity of ^' /a' is one of the conditions for the stability. Another 
condition is the positivity of the hi parameters defined in eq. (j22I)- We can 
easily see how 62 changes with a. From definition ()22|) and equation (0) with 
the V^'-term neglected, one obtains 



The exact behavior of the r.h.s. of this equation depends on the details of a 
model, but one feature is rather model independent. We have shown that the 
warp factor evolves faster for bigger values of a. This means that typically 
62 increases with increasing a. 

The main consequences of adding the Gauss-Bonnet term to the action 
are discussed in the two previous paragraphs. Such a term with a positive 
coefficient a causes the inter-brane distance to decrease, and the the radion 
mass squared to increase (improving the stability of the brane positions). The 
results of adding the GB term with a negative coefficient are more model- 
dependent, but in general the stability seems to be worse (the radion becomes 
lighter, and can eventually become tachyonic). 

Let us now discuss the inter-brane distance in more detail. Of course the 
second brane can be placed at 1/2 smaller than the maximal value discussed 
after eq. ()56|). A very important point is that the boundary conditions © 
and ()10|) must be fulfilled at that second brane. The values of the l.h.s. 
of both the equations are, for a given background solution, known at each 
y. For a fixed form of the potential U2-, also the r.h.s. of these equations 
are known for each y. For a general value of y none of those two boundary 
conditions is fulfilled. The idea of the Goldberger-Wise mechanism j3] is 
that there is a special value of y for which both equations and (fTUj) 
are fulfilled, and the inter-brane distance is dynamically determined to be 
y2 — Hi- Of course, in general this can not be achieved. The reason is very 
simple: there are two conditions but only one parameter to be adjusted, 
namely y2- So we need a mechanism which adjusts the parameters of the 
Lagrangian in such a way that one combination of the boundary conditions 
is fulfilled "automatically". This is just another form of the cosmological 
constant problem. Unfortunately, the result given by the Goldberger-Wise 
mechanism depends very strongly on the details of the unknown solution to 
that problem. In the next section some results of numerical calculations are 
presented to illustrate our analytical analysis. 
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7 Numerical examples 



We consider a simple model similar to that discussed in [0]. The bulk and 
brane potentials have the following forms 

V{ct>) = A + ^M202 , Ui{ct>) = \ + - . (59) 

We chose the parameters in V and Ui as follows: A = —15, M = 1, Ai = 10, 
/ii = 100, t>i = 1 (in K = 1 units). For 5 values of a, we solve the boundary 
condition at ?/i = numerically, and integrate the background equations 
of motion. Two stability conditions, 6i > and (t)'{y) 7^ 0, are fulfilled in 
all these cases. The parameter ^' /a' as a function of y is shown in fig. 1. 
The results are consistent with our qualitative discussion: away from the 
first brane, ^' /a' is bigger for bigger values of a. This effect increases with 
y. For positive a, the ratio ^' /a' becomes very large when we approach 
the singularity of the background solution. For the Einstein-Hilbert theory 
(a = 0), this ratio decreases but stays positive up to the singularity. For 
negative a, it decreases faster and becomes negative before the corresponding 
singularity. 




0.0 0.1 0.2 0.3 0.4 



y 

Figure 1: Ratio E,' /a' as a function of y for different values of a: 0.01 (a); 
0.005 (b); -0.005 (c); -0.01 (d). The curve without a label corresponds to a 
model without the Gauss-Bonnet term (a = 0). 

Then, for each value of y, we solve the boundary conditions numerically as 
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if the second brane were positioned at that y. There are two such boundary 
conditions and three parameters in U2 so one of the parameters can be fixed. 
First, we fix A2 = — Ai in order to compare our results with those of jHj- The 
obtained values of 62 are plotted in fig. 2. The model can be stable only for 
positive 62, so the second brane can be placed only at the values of y for 
which 62 > 0. 




0.0 0.1 0.2 0.3 0.4 



y 

Figure 2: Parameter 62 defined in eq. fj22j) as a function of the inter-brane 
distance for A2 = — Ai. Labels are as in fig. ^ Curves for a > end at the 
background solutions singularities, curves for a < end at values of y for 
which ^' /a' becomes negative. 

We have discussed the values of ^' /a' and 62- In section 4, it has been 
shown how they are related to the radion mass. To illustrate this relation, 
we have calculated the radion mass numerically for some cases. For the 
parameters given above and for y2 — yi = 0.15, we have found that ttiq equals 
-0.12, 0.07, 0.40, 1.01, 2.43 for 5 values of a: -0.01, -0.005, 0, 0.005, 0.01, 
respectively. Comparing this with figs. ^ and |21 one can see that indeed ttiq 
grows with the values of C,'/a' and 62- A negative value of rriQ for a = —0.01 
reflects the fact that in this case 62 < 0. 

We confirm the result of [S] that in this particular set-up there is a min- 
imal inter-brane distance below which the model is unstable. With GB 
interactions, this minimal distance decreases (increases) for positive (nega- 
tive) a. But there is also an upper bound on the inter-brane distance. For 
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a > 0, it is determined by the singularity of the background solution^ while 
for a < it is given by value of y for which ^'/a' becomes negative. 

One should ask the question: is the above result typical for the considered 
model? The answer is: no, it is not. To illustrate this, we show in figs. 3 
and 4 the y-dependence of the parameter 62 for models where the value of 
A2 is changed by 10%. Such small changes in one of the parameters have an 
important influence on the possible inter-brane distance. There is no lower 
bound on y2 ~yi for A2 = — 0.9Ai or for A2 = — l.lAi. But in the second case 
not all values of y2 — yi between and the maximal value are allowed. There 
is a window of forbidden values of y for which the model becomes unstable 
while being stable outside. 

In fact, the functional behavior of 62 shown in figs. El and |3] is much more 
typical than that in fig. |21 For general values of A2 there is no lower bound 
on y2 — yi- For A2 ~ — Ai, there are two ranges of y2 — yi for which the model 
can be stable. The upper range shrinks when we decrease A2, and eventually 
disappears for large enough IA2I. The lower bound on y2 — yi exists only 
when A2 ~ — Ai up to a few percent. The common feature of all the models 
is the existence of an upper bound on the inter-brane distance above which 
the radion can not be stabilized. 

The existence of such an upper bound can cause problems for model build- 
ing. The most appealing rationale for considering the 5-dimensional brane 
models comes for the Hofava-Witten construction motivated by M-theory 
p. One of the main features of those models is a relatively big length of 
the 5-th dimension necessary to obtain a large enough 4-dimensional Planck 
scale. Hence one has to check whether an upper bound on the inter-brane 
distance in a given model is compatible with phenomenological constraints. 
This is especially important after including GB interactions, which typically 
cause the maximal inter-brane distance to decrease. 

We would like also to stress (again) that radion stabilization is quite 
sensitive to the unknown mechanism of solving the cosmological constant 
problem. In the examples presented in this work we fixed A2 and used two 
other parameters in U2 to fix the brane position and to set the cosmological 
constant to zero. Of course, we do not know whether these two particular 
parameters have anything to do with the cosmological constant problem. 
Hence one can consider situations where parameters other that 112 and f 2 are 
adjustable. We have checked that, for example, fixing V2 gives three classes 
of solutions for 62 similar to those shown in figs. |21 to 01 However, details of 
a given solution depend on the actual value chosen for t>2. 

^Such an upper bound exists also in the absence of the Gauss-Bonnet term, which was 
not observed in ref. [HI 
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8 Conclusions 



We have found the equations of motion for scalar perturbations in the pres- 
ence of the Gauss-Bonnet interactions for an inflating background. The 
expression for the mass of the lightest of such perturbations, which can be 
identified as the radion, has a form very similar to that for the standard 
Einstein-Hilbert gravity. However, the presence of the Gauss-Bonnet term 
with a positive coefficient typically increases the actual radion mass. This 
can help to stabilize the radion in an inflating background. 

The necessary and sufficient conditions for the radion stabilization have 
been found for the case of a vanishing Hubble constant. One additional con- 
dition appears when GB interactions are taken into account. The form of 
the other conditions does not change after adding those interactions. How- 
ever, the parameters of the model for which those conditions are fulfilled 
can change substantially, because the character of the background solutions 
changes when the GB term is present in the action. 

There are two classes of solutions in the presence of GB interactions. In 
the limit where the GB coefficient goes to zero, solutions from one class just 
converge to the solutions in the Einstein-Hilbert theory, while solutions from 
the second class do not. It has been shown that the solutions from the second 
class are strongly disfavored by a stabihty analysis. 

Even without inflation, numerical calculations are in general necessary to 
check whether the branes can be stabilized for a given model. The results for 
the radion stabilization are quite sensitive to all (mainly unknown) details of 
models. They depend, for example, on assumptions made about the mecha- 
nism leading to an acceptably small cosmological constant. This is true for 
models both with and without the Gauss-Bonnet interactions. 

However, some typical features of the solutions can be found using ana- 
lytical arguments. The most important of them is the existence of an upper 
bound on the inter-brane distance for each given model. General changes 
caused by the GB term can also be identifled. Typically, the values of some 
parameters appearing in the formula for the radion mass change with a in 
such a way that this mass grows with growing a. Hence (with all other pa- 
rameters fixed) stability increases if the GB term with a positive coefficient is 
present in the action. At the same time, the maximal value of the inter-brane 
distance becomes smaller. 
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